Phase diagram of the half— filled Hubbard chain with next— nearest— neighbor hopping 
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We investigate the ground-state phase diagram of the half-filled one-dimensional Hubbard model 
with next -nearest-neighbor hopping using the Density-Matrix Renormalization Group technique as 
well as an unrestricted Hartree-Fock approximation. We find commensurate and incommensurate 
disordered magnetic insulating phases and a spin-gapped metallic phase in addition to the one- 
dimensional Heisenberg phase. At large on-site Coulomb repulsion U, we make contact with the 
phase diagram of the frustrated Heisenberg chain, which has spin-gapped phases for sufficiently 
large frustration. For weak U, sufficiently large next-nearest-neighbor hopping t<z leads to a band 
structure with four Fermi points rather than two, producing a spin-gapped metallic phase. As 
U is increased in this regime, the system undergoes a Mott-Hubbard transition to a frustrated 
antiferromagnetic insulator. 
The one-dimensional Hubbard model is the prototyp- 
ical model for strongly interacting electrons in one di- 
mension. For repulsive interaction, its low-energy, long- 
distance physics is well-described by the Luttinger liquid 
picture, in which the fundamental excitations are gap- 
less bosonic spin and charge modes, and the correlation 
functions .exhibit critical behavior with non-universal 
exponentsu. At half filling, Umklapp processes lead to a 
gap in the charge excitation spectrum and thus insulat- 
ing behavior for any finite value of the on-site Coulomb 
interaction, U. The spin excitations behave as in the 
strong-coupling, Heisenberg limit, i.e. are gapless with 
linear dispersion. 

The introduction of a next-nearest-neighbor hopping 
can change this picture dramatically. In strong cou- 
pling, the additional hopping leads to a frustrating 
next-nearest-neighbor Heisenberg interaction so that the 
model maps to the frustrated Heisenberg chain. At weak- 
coupling, the effect of ti is to change the band structure, 
and, in particular, the number of Fermi points. In this 
paper, we shall explore the interplay between the frustra- 
tion at strong coupling and the changed band structure 
at weak coupling. As we shall see, the resulting phase 
diagram contains a number of highly interesting phases: 
a spin-gapped metallic phase, commensurate and incom- 
mensurate disordered magnetic insulating phases, as well 
as the one-dimensional Heisenberg insulator. 
We study the Hamiltonian 



FIG. 1. The ti-t 2 Hubbard chain. 

Here the summation goes over L sites and spin a, and 
we will always take U positive and set ti = 1. Since 
the sign of ti is irrelevant at half filling due to particle- 
hole symmetry, we only consider ti > in the following. 
For U = and periodic boundary conditions, H can be 
diagonalized via a Fourier transform, yielding 
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with k an integer multiple of -f- and 
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-2ii cos k — 2ti cos 2k . 
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An interesting feature of this band structure is that there 
is a nontrivial transition as a function of t% even at U = 0. 
For ti < 0.5, the noninteracting band has two Fermi 
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points and for ti > 0.5, it has four Fermi points (±k 

(2) 

and ±k F ). This is important in a weak-coupling picture 
because the Fermi points are separated by the Umklapp 
vector q = it only for t% < 0.5. 

For large U, Eq. ([!]) can be expanded perturbatively in 
1/U, leading to the one-dimensional frustrated Heisen- 
berg Hamiltonian 
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H — (J%SiSi+i + JiSiSi. 
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where t\ is the nearest-neighbor and ti the next-nearest 
neighbor hopping and U is the on-site Coulomb repulsion. 
It is useful to visualize the geometry as a zigzag structure 
as depicted in Fig. [I]. 



with Ji = -jj- and Ji = -jf-. This model has been ex- 
tensively studied using a number of different methodso. 
In particular, there is an exact solution at JilJ\ = 0.5 
(tijti = 1/a/2) due to Majumdar and GhosbE. At this 
point, the ground state is a simple dimer configuration 
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and the spectrum has a spin gapa For J2/J1 < J c ~ 
0.241167, studies combining numerical and field theory 
calculationsO have shown that the spin excitation spec- 
trum is gapless. In the strong-coupling expansion, this 
parameter value maps to the point at which the Fermi 
surface jumps from two points to four in the t\-t2 Hub- 
bard chain. When J c < J2/J1 < 0.5, dimerization cor- 
relations are present ancLfor J2/J1 > 0.5, i.e. above the 
Majumdar-Ghosh pointu, incommensurate spiral spin 
correlations appear. 
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FIG. 2. The unrestricted Hartree-Fock phase diagram. 
The insulating phases are labelled by the wave vector q of 
the magnetic ordering for the phases in which the order pa- 
rameter m > 0. 

Mean-field theory can be a useful tool to provide quali- 
tative information about the ground-state phase diagram 
even though it tends to overemphasize ordered phases 
and cannot yield the correct critical behavior. Here we 
perform mean-field calculations starting from the unre- 
stricted Hartree-Fock Hamiltonian 
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where H is the non-interacting Hamiltonian, ng is the 
total electron density on a site, and se = s , c\ s a s ^ s iC£ S > 

with (Tj * z ' the Pauli matrices. The mean fields pe and 
me are determined by the self-consistent equations pi = 
(ne) and = (s^). We postulate a uniform density 
pe = p, so that J2ePe n £ = P-^ is constant, and a spiral 
arrangement of the magnetic moment 



mg = m (cos q£, sin q£, 0) . 
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The ground-state energy Eq is then a function of the 
variational parameters q and m and the mean-field gap 
is A = ^p. In order to obtain q and m, we minimize 
Eo(q,m) numerically for large systems. When m = 0, 
the system is metallic, but if m is finite, the half-filled 
system is a magnetically ordered insulator with principal 
wave vector q. 



We obtain the mean-field phase diagram shown in Fig. 
0. For t% < 0.5, the system is an antiferromagnet (i.e. 
q = 7r) for all U > 0. This phase is the same as that ob- 
tained for £2 = 0: The Fermi points are separated by the 
wave vector q = n, leading to antiferromagnctic ordering. 
For t2 > 0.5, this commensurate separation of the Fermi 
points is absent, and the system is a paramagnetic metal 
for small U. 

As U is increased, the system undergoes a transition 
to an insulating phase with magnetic ordering at wave 
vector q = k^ — kp = n/2. For still larger U, there is 
an incommensurate phase with ordering at wave vector 
q* , where q* goes continuously from it at smaller £2 to 7r/2 
as t2 — * 00. As we shall see in the following, the phase 
boundaries and wave vectors found in this mean-field 
phase diagram are qualitatively similar to those found for 
the fully interacting system using the DMRG, although 
the nature of the phases themselves is different. 
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FIG. 3. DMRG phase diagram at half filling. The symbol 
□ indicates a C0S1 phase, • a C0S0 phase and * a C1S0 phase. 
The approximate phase boundaries are marked as a guide to 
the eye. 

Using the DMRgI, we can obtain the ground-state en- 
ergy and wave function on a finite lattice numerically to 
very high accuracy. Here we perform calculations keeping 
up to 800 states on lattices of up to 64 sites so that the 
maximum weight of the discarded density matrix eigen- 
values is 10~ 6 . In order to determine the phase diagram, 
we first calculate the charge and spin gaps, defined as 

A p = ~[E (N + 2,0)+E (N-2,0)-2E (N,0)} (7) 
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E (N,1) - E o (N,0) , 
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where Eq(N,S) is the ground-state energy for N parti- 
cles and spin 5*. We calculate the excitation gaps for 
systems with different size L and then extrapolate to 
L — > 00 using a quadratic polynomial in We have 
calculated the spin gap as a function of U in the large 
U limit (e.g. at U — 100) and find very good quantita- 
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tive agreement with numerical results for the frustrated 
Heisenberg chain from Ref.a. 

In Fig. pL we display the phase diagram in the U— i 2 
plane determined using the L — > oo spin and charge gaps 
computed with the DMRG. We label the phases using 
the notation of Ref. [?], in which CnSm represents a phase 
with n gapless charge modes and m gapless spin modes. 
For ti < 0.5, the system is in the one-dimensional Heisen- 
berg phase found for t 2 = 0, a C0S1 phase. As in the 
mean-field phase diagram, this phase extends to the en- 
tire region in which there are two Fermi points. For 
t 2 > 0.5, there are two different regions both with gapped 
spin modes. For weak U, the system is metallic (C1S0), 
and for U > U c of the order of the bandwidth, the sys- 
tem is insulating (C0S0). In this region, where the sys- 
tem has four Fermi points, a weak-coupling renormal- 
ization group calculation predicts a metallic spin liquid 
i.e. C1S0, phase,H. Such a spin-gapped phase has been 
found numerically for this model away from half fillingta. 
At half filling, there are no relevant Umklapp processes 
which could drive the system to an insulating phase when 
there are four Fermi points!! 
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FIG. 4. Charge gap as a function of U for t% 
filling. The dotted line is linear in U. 
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In Fig. ||, we show the charge gap as a function of U for 
t 2 = 1. There is a clearly defined metal-insulator tran- 
sition at U c — 3.2 w W/2, where W is the bandwidth 
of the non-interacting case (W = 6.25 for t 2 = 1). For 
U ^> U c , the gap grows linearly as U— LL as one would ex- 
pect above a Mott-Hubbard transition!!], but the curve is 
rounded near the transition point. The metal-insulator 
transition as a function of U therefore appears to be con- 
tinuous. 
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FIG. 5. (a) Fourier transform, S(q), of the spin-spin corre- 
lation function for the half-filled system with L — 40, U = 100 
and ti — 0.7 . . . 1.5 in steps of 0.2 (solid lines from bottom to 
top at q = tt/2); the dashed line is the result for — 0. (b) 
The wavevector q* obtained from the position of the peak in 
S(q) as a function of ti for various values of U. 

We now investigate the nature of the magnetic order- 
ing in the different phases by examining the spin-spin 
correlation function 



Sa. y (r) 



S t S £+r) 



(9) 



where (sj) are the spin raising (lowering) operators 
corresponding to se, and we average over a number of 
^-values (typically six) to reduce oscillations present be- 
cause of the open boundaries. We then perform a con- 
tinuous Fourier transform to obtain the static structure 
factor 
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The resulting function is plotted in Fig. ||(a) for various 
i 2 ranging from 0.7 to 1.5. It is clear that the peak shifts 
continuously from q = it at i 2 = 0.7 to q = tt/2 at 
t 2 = 1.5. As t% increases, the peak at q = tt/2 becomes 
sharper, resembling more closely the q = n peak of t 2 = 
structure factor, also shown for reference. This is because 
the correlation length diverges as < 2 becomes largeo; it is 
the limit of two uncoupled Heisenberg chains. 
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The position of the peak, q* , is plotted as a function 
of t 2 in Fig. |(b) for different U values. For U > 5, 
q* = 7r for t2 < t\ ~ 0.7. At large U, one expects q* to 
deviate from tt oiiIm for <2 above the Majumdar-Ghosh 
point, tf G = I/V2A. As can be seen, t\ m tf G down to 
small values of U. Note that t* 2 does not coincide with the 
opening of the spin gap in Fig. ^ there is an intermediate 
region of a spin-gapped dimer phase with q* = it. AsU is 
reduced, the region of incommensurate spiral spin order 
becomes narrower. This is expected because there is no 
incommensurate spiral spin order at U = 0. 
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FIG. 6. Pairing correlation function on a log-log scale for 
L = 64, t 2 = 0.78 and U = 0, 1. 

In the metallic spin-liquid (C1S0) region, the long- 
distance behavior is characterized by the correlation 
functions associated with the charge degrees o£ freedom. 
A bosonization treatment valid for weak [70 predicts 
that two competing correlation functions have the slowest 
asymptotic decay: the dimer-wave correlation function 
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and the singlet superconducting correlation function 
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If 9 > 1, the pairing is dominant, while for 9 < 1 the 
dimer wave is dominant. In addition, at half filling high- 
order Umklapp process are relevant for 9 < 1 and the 
system becomes insulating!! We then expect that d-wave 
pairing correlation are dominant for U < U e , where the 
charge gap opens at U c . We compare the pairing function 
for the weakly interacting and the non-interacting case to 
see if there is an enhancement. At half filling kp^ — fc£P = 
7r/2, so the relevant pair operators in real space involve 
pairs on next-nearest-neighbor sites. We have carefully 
evaluated the pairing correlation function 
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where Aj = (c^Ci + 2[ — CiiCi+2f) / \/2 and we average over 
a number of i-values to reduce oscillations present be- 
cause of the open boundaries. In Fig. ^, we plot the re- 
sults for U = and {7 = 1, and we clearly see that there 
is no enhancement when we turn U on. This_is in contra- 
diction with previous work by Kuroki et al.llil who claim 
that the superconducting correlations are dominant for 
small U based on projector QMC. We have also evalu- 
ated the dimer-wave correlation function which shows no 
enhancement for U < U c , but the long range behavior 
clearly changes from r~ 2 to when U > U c . 

In conclusion, we have explored the rich ground-state 
phase diagram of the half-filled t\— t-x Hubbard chain 
using the Density Matrix Renormalization Group. For 
t2 < 0.5 the model is insulating for all values of U due to 
Umklapp processes and there is no spin gap. For t% > 0.5 
the system has a spin gap for all U but is metallic (i.e. 
the charge gap vanishes) for sufficiently small U. A con- 
tinuous metal to insulator transition occurs when U is of 
the order of the bandwidth. At large U, we find quantita- 
tive agreement with the phase diagram of the frustrated 
Heisenberg chain. While the overall phase diagram is 
in agreement with weak-coupling renormalization group 
calculationsa, the pairing correlations do not behave as 
predicted in the spin-gapped metallic phase. The ex- 
act nature of the metallic phase is difficult to determine; 
more work must be done in this direction. 
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